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THE REDUCTION OF TOPSIDE IONOGRAMS TO 
ELECTRON DENSITY PROFILES 

by 

John E. Jackson 
Laboratory for Space Sciences 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

The basic principles of the ionosonde technique 
are reviewed with special emphasis given to the satel- 
lite version of this technique (topside soundings) . 

It is shown that the lamination method (developed ini- 
tially for the analysis of ionograms obtained with 
ground-based ionosondes) can be used to calculate the 
electron density N vs. true height h in the topside 
ionosphere provided a number of refinements are intro- 
duced. These refinements are required because the 
extraordinary wave data has to be used to obtain N(h) 
profiles from topside ionograms, and because these data 
are very sensitive to the magnetic field variation 
over the great altitude range of topside soundings 
(typically one order of magnitude greater than for 
ground-based soundings) „ A complete procedure for 
analyzing topside ionograms, based upon a parabolic- 
in-log N lamination method is presented which incor- 
porates the actual values of the earth magnetic field 
at all heights, a change in variable which renders the 
group path integrand finite at the reflection point, 


i 



(and varying sufficiently slowly elsewhere to be cal- 
culated very accurately with a 3-point Gaussian inte- 
gration technique), and iteration until successive 
calculations agree to within 0.01 km. A procedure was 
also developed to insure and accelerate the convergence 
of the iteration process. The importance of each re- 
finement is discussed and illustrated quantitatively. 

It is shown that this new N(h) reduction technique can, 
in theory, yield topside electron density profiles with 
an altitude accuracy of a few kilometers over a 3000 to 
400 km altitude range. )( This theoretical accuracy can be 
achieved using typically ten to twenty data points per iono- 
gram. In actual practice the errors in the N(h) profiles 
arise primarily from the limited resolution of the iono- 
grams, calibration inaccuracies, and deviations from 
vertical propagation. The actual accuracy of topside 
N(h) profiles is discussed in a companion paper (Jackson, 
same I.E.E.E. issue). 
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THE REDUCTION OF TOPSIDE IONOGRAMS TO 
ELECTRON DENSITY PROFILES 

by 

J. E. JACKSON 


I, INTRODUCTION 

Most of our knowledge concerning the distribution of electrons 
in the ionosphere as a function of position and time has come from 
ionospheric sounders (ionosondes) located either on the ground or 
in satellites. Ground-based soundings provide the electron density 
N as a function of altitude h below the altitude of maximum 
electron density hmaxF2 , (typically at 300 km) and topside 
soundings provide the N(h) data from the topside sounder altitude 
down to hmaxF2, Other experimental configurations which have been 
used include sounders in aircraft and in rockets. Thus sounders 
can be either mobile or geostationary and be located either above 
or below hmaxF2. 

Ionosondes operate on principles similar to those of radar. 
There are, however, important differences. The ionosonde does 
not receive echoes from hard targets as is usually the case with 
radar. The ionosonde target is a plasma, whose reflection pro- 
perties are highly frequency-dependent. The reflection conditions, 
which will be discussed later, are such that 1) the ionosonde must 
utilize frequencies lower than those normally used in radar and 
2) a large number of sounding frequencies are required in order 
to investigate the plasma and obtain its electron density distri- 
bution. This last requirement is best met with the swept-frequency 
system which is commonly used in ground-based ionosondes, and which 
was selected for the Alouette I and Alouette II sounders. Single 
or multiple fixed-frequency sounders have also been utilized for 
special applications. For example, the Explorer XX sounder operated 
at six fixed frequencies (references to this issue to be provided 
by editors) . 

In conventional radar applications the round trip time for a 
reflected signal is readily converted to target distance, since the 
radar signal travels essentially at the velocity of light in vacuo. 
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A similar round-trip-time to distance conversion is required in 
the reduction of the ionosonde data, but in this case the problem 
is much more difficult . The conversion process must take into 
consideration the fact that the signal velocity varies continuously 
while the props gation path is within the ionosphere . In the 
ionosphere the signal velocity is less than the velocity of light 
in vacuo resulting in a retardation of the signal. The calculation 
of this retardation, so that real distances to the point of reflection 
can be determined, is a subject which has received the attention of 
ionospheric workers since the earliest days of radio sounding 
experiments. The techniques used to perform these calculations 
are based upon the magneto-ionic theory (Budden 19G1, Kelso 1964, 
Ratcliffe 1959). A single reduction technique, however, cannot be 
used for all the possible experimental configurations, because each 
presents its own special problems. Thus techniques developed during 
the 1940 and 1950 decades for the reduction of ground-based ionograms 
were not directly applicable to the analysis of topside ionograms. 

For a number of reasons, including practical considerations 
such as the need for minimizing computation time, a large number 
of techniques (Radio Science, Oct. 1967) were developed for the 
reduction of ground-based soundings. The availability of fast 
computers has led to the adoption of the lamination technique 
(discussed later) for most of the N(h) work based upon ground-based 
soundings. In view of its flexibility and accuracy, the lamination 
technique was also selected for the analysis of topside ionograms. 
Extensive modifications and refinements had to be incorporated in 
this technique in order to meet the special requirements of topside 
N(h) analysis. Much of the present paper is devoted to a discussion 
of these topics. 

The ionogram reduction techniques in common use all assume 
that the soundings are vertical. This assumption, in general, does 
not introduce serious errors in the reduction of ground-based 
soundings. For topside soundings, however, the assumption can 
lead to large errors in N(h) results. The magnitude of the error 
due to non-vertical propagation tends to increase with satellite 
altitude. On Alouette II, for example, the ionograms taken from 
altitudes above 2000 km are frequently the results of reflections 



obtained under conditions of oblique propagation. Oblique 
propagation is much less of a problem on Alouette I ionograms 
which corresponds to soundings taken from a nearly constant 
altitude of 1000 km, or on Alouette II ionograms when the Alouette II 
altitude is less than 1500 km. Proper care must therefore be 
exercised in the selection of topside ionograms for N(h) analysis, 
and when this is done the errors in the derived profiles are usually 
less than 10 percent. This error estimate is based upon typical 
results obtained when profiles derived from topside ionograms were 
compared with profiles obtained simultaneously by other techniques, 
such as rocket and incoherent backscatter measurements (see paper 
in this issue "Comparisons between topside and ground-based 
Soundings" by J. E. Jackson). 

Using procedures more elaborate than those employed for 
vertical soundings, it is often possible to derive reliable N(h) 
profiles from ionograms obtained under conditions of oblique 
propagation. For example, consecutive high altitude ionograms 
often show clear traces corresponding to propagation along magnetic 
field lines; these traces can be used to compute a series of field 
aligned profiles from which the N(h) profiles can be derived (Colin, 
Dufour and Willoughby, this issue). In this case the initial field- 
aligned profiles are calculated using the same basic methods des- 
cribed here. With more complex reduction techniques (Lockwood, 
this issue), vertical N(h) profiles can in many cases be obtained 
from more general (non field-aligned) oblique soundings. These 
methods, however, are not presently in general use and they are 
not the subject of this paper. 

Although some of the topside ionograms are unsuitable for 
N(h) analysis, they often provide the basis for investigating 
various ionospheric phenomena such as partial field-aligned ducts, 
spread-F, and other irregularities (reference to appropriate papers 
in special issue) . Another important use of ionograms has been for 
the study of plasma resonances (reference to special issue) . The 
topside sounder is therefore a very versatile tool for ionospheric 
studies, the N(h) analysis discussed in this paper being only one 
of the many applications of topside soundings. 
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II. DISCUSSION OF p’-f DATA ON TOPSIDE IONOGRAMS 


The ionograms which will be discussed are those obtained with 
topside sounders which sweep in frequency from about 0.1 to 20 MHz 
emitting short pulses at 20 to 60 pulses per second. The received 
echoes are conventionally displayed in the "ionogram" format 
(Franklin, et al; Hagg; this issue) in which the time delay between 
pulses is displayed as a vertical axis, calibrated in terms of 
apparent range p ? . This calibration assumes that the v> a ve travels 
at the velocity, c, of light in vacuo. In general p’ is larger than 
the actual distance to the echoing region because the sounding 
signal is retarded by the ionosphere. By calculating the extent 


to which each sounding frequency is retarded, the N(h) profile 
can be computed. The group velocity V G of the sounding signal 
varies as a function of frequency f , electron density, terrestrial 
magnetic field B and the angle cp between B and the direction of 
propagation. It is the complexity of the function which has made 
it necessary to develop special techniques for converting ionograms 
to N(h) profiles. 

By definition, the apparent range at a frequency fi is one 
half the measured round-trip time (Ati) multiplied by the velocity 
of light in vacuo. Consequently: 

P'(fl) = (Atl ) QX 


The round-trip time is given by: 

Ati — 2 


Pi 
P « 


dp 


( 2 ) 


where the integration is along the path taken by the wave, p Q is 
the position of the sounder and Pi corresponds to a reflection point 
for the frequency fi . Combining equations (1) and (2) gives: 


P'(fi) 


f>Pi c 

1 (y^) dp 

p o G 


(3) 


- 4 - 



In standard magneto-ionic nomenclature the quantity c/Vq is called 
the group refractive index n'. Equation (3) then becomes: 


P* 



n ’dp 


(4) 


* 

The behavior of n’ is usually described in terms of a parameter X 

which is proportional to N/f 3 . Although this conventional approach 

is used later in the text . it is more convenient here to use a V_ - 

G 

vs - f diagram. A V^-f diagram is shown in Figure 1 for one set 
of values of N, B and cp. In spite of broad variations, certain 
features are common to all V^-f diagrams. F-rst , there are always 
three group- velocity curves as indicated on Figure 1. These curves 
are conventionally called Z, ordinary (0) and extraordinary (X) . 

The "ordinary" mode is the only one which would be present in the 
absence of a magnetic field. Second, there are two values of Vq 
over most of frequency range for which propagation is possible (V^ > 
They correspond to the cwo independent modes of propagation that 
are possible when B / 0, the ordinary and the extraordinary modes. 
These two modes are present on the ionograms because a plane 
polarized wave is transmitted by the sounder and its two circularly 
polarized components propagate independently at their respective 
group velocities. The reflection points for the Z, (0) and (X) 
components are labelled fz, fN, and fx. The plasma frequency fN is 
a function of N only, whereas fz and fx depend upon both N and B. 

The apparent ranges for the Z, (0) and (X) traces are related to 
the true range at any given point by an integral such as the one 
shown in Equation (3). The function appropriate for the parti- 
cular mode must of course be used in the integral. The point on 
the Z curve labelled fzl (fz, infinity) is called the Z-mode cut 
off and is the maximum frequency at which the Z wave can propagate 
for the given values of N, B, and v. 
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The curves of Fig. 1 were calculated for conditions existing 
at the satellite when the topside ionogram of Fig. 2 was obtained. 
Thus, Fig. 1 indicates the frequencies at which fz, fN, and fx 
would be seen on the ionogram for conditions near the satellite. 

The letter S is often adoed to fx and fz as done in Fig. 2 when 
fx and fz refer to conditions near the satellite. On topside 
ionograms the Z- mode cut-off fzl occurs at the satellite, resulting 
in very large Z-wave delays as the sounding frequency approaches 
fzl, and consequently the fzl phenomenon is exhibited as a maximum 
range effect for the Z-wave. The complete Z, (0) and (X) traces 
are of course influenced by the gradual changes in ambient conditions 
along the propagation path, i.e., by the continuous changes in 

v G (N, b, cp » f) fro.a the satellite altitude down to the various re- 
flection points. There are other features on the ionogram, such 
as the resonances occurring at fH (gyrofrequency) and at fT (upper 
hybrid frequency). A strong resonance is also present on this 
particular ionogram at the second harmonic of the gyrofrequency (2fH) . 
Plasma resonances are discussed in other papers of this issue 
(appropriate references from editors). Also visible on Fig. 2 are 
echoes due to the side-band responses of the transmitted pulse 

(See Hagg & Hewens, Warnock, this issue). There is one response 
below the main spectrum and two responses above. These si-de-bands 
are very clearly seen on the Z’ trace of Fig. 2. This Z T trace is 
an additional Z trace originating from the plasma resonance 
(Calvert, 1966). 

In principle, electron density profiles can be calculated 
from either the Z, (0) or (X) traces. Since the three traces are 
affected quite differently by the ionospheric conditions, the 
analysis techniques can be checked by comparing results obtained 
from each of the three traces. The results of such a comparison 
using the data of Fig 2 are shown in Fig. 3. However, as seen in 
Fig. 1, the Z-wave is normally cut-off at the satellite and does 
not penetrate down to the maximum of the F-2 region. Since ex- 
ceptions to this rule are extremely rare (Jackson, 1967, p. 21), 
the Z-trace (Colin, this issue) is seldom used for N(h) analysis. 
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This difficulty, however, is not present on either the (0) or the 
(X) traces, and either trace can be used to derive an N(h) profile. 

An analysis based upon the ordinary trace is somewhat simpler, 
because the reflection density is independent of the magnetic field. 

A more important consideration, however, is the relative quality 
of the (0) and (X) data. 

In the lower ionosphere (D and E regions) , collisions between 
electrons and neutral particles cause low frequency waves to be 
attenuated. The absorption is more severe for the (X) mode than 
for the (0) mode. Consequently on ground-based ionograms the (0) 
trace is usually more complete and therefore it is the one used in 
N(h) analysis. The situation is different on topside ionograms. 
Attenuation due to collisions is negligible; the controlling factors 
are now the sounder antennas and the fact that reflection at a 
given density N R occurs at a higher frequency for the (X) mode 
than for the (0) mode. For example, if N R = 1240 el/cc and 
fH = 1.0 MHz, the (0) ray reflects at f = 0.316 MHz and the (X) 
ray reflects at f = 1.09 MHz. Although the antennas used in a top- 
side sounder satellite are physically very long, they are electrically 
short (and hence difficult to match to the transmitter) at the low- 
frequency end of the sweep. Based upon this consideration alone, 
transmission at 0.316 MHz would be considerably weaker than at 
1.09 MHz. The situation is made even worse because the antennas 
are immersed in the ionosphere. This situation causes a change 
in antenna impedance, which is particularly severe near the plasma 
frequency. The net result is that the low-frequency end of the (0) 
trace is often missing on topside ionograms, and consequently the 
(X) trace is normally used for the N(h) analysis. 

The preceding qualitative discussion omitted most of the usual 
formulas. The quantities of interest and the corresponding formulas 
are listed below. The basic ionospheric parameters which affect 
radio soundings are: 

3 

N = ionospheric electron number density (electrons/cm ) 

B = induction (gauss) of terrestrial magnetic field 
cp = angle between B and direction of propagation 
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The plasma frequency fN (ordinary ray reflection point) and the 
gyrofrequency fH depend only upon one ionospheric parameter as 
indicated below: 


fN (MHz)= (8 - 98xl0 ' 3 > /N 

(5) 

fH CMHz)=< 2 - 8)B 

(6) 


The upper hybrid frequency fT, the extraordinary wave reflection 
point fx and the Z wave reflection point fz depend upon both N 
and B as follows: 


fT= \ fN 2 +f H 2 

/ 

(7) 

, fH . V4fN 2 +fH 2 
lx T + 2 

(8) 

f z-f x— fH 

(9) 


From Equation 8 and 9, 

fN 2 =fx(fx-fH)=fxfz = (N/12,400 if frequencies are in MHz) (10) 


The high frequency cut-off (fzl) of the Z-wave depends upon N, 
B and cp as indicated below: 


fzl 


fT 2 + V fT 4 - 4f N 2 f H 2 COS 2 cd 


( 11 ) 


Finally, the group velocity V^, (or equivalently the group refractive 
index n^c/V^) and consequently also the group path integral 
Pi = f^n’dp depend upon not only N, B and cp, but also upon the 
frequenSy of the sounding signals. 

For the special case in which reflections are vertical, p = h, 
and cp = 90°-9, where 0 is the magnetic dip. Since vertical pro- 
pagation is assumed in the subsequent analysis, the group path 
integral will be written in terms of real heights as follows: 
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h R (f) 

p'(f) = r n’[N(h), f ,B (h) , 6 (h) ] dh (12) 

h 0 

where : 

h = actual height 

h^ = height of sounder (satellite altitude for a topside 
sounder) 

h R = height at which reflection occurs for frequency f 

n' = group index of refraction 

The basic problem which must be solved in the N(h) analysis 
is the inversion of the group path integral, i.e., the conversion 
of the p' (f) function to the N(h) function. The remainder of this 
report discusses this conversion for the case of vertical propagation. 
It should be noted that Equation (12) does not give the function n' 
in explicit form. The formula for n’ is quite complicated, (see 
appendix A) and is not needed for the present discussion. Tables 
giving values of n' under a wide range of conditions have been 
published (Becker, 1960). 

Ill . OUTLINE OF THE LAMINATION CONCEPT 

For a given geographic location and a given N(h) distribution, it 
is a relatively straight forward matter to evaluate the integral 
shown in Equation (12). For a given frequency f, the density N R at 
the reflection point is known from Equation (10) , and the integration 
limit h R follows immediately from the known N(h) function. The 
magnetic field parameters are known as a function of altitude. 

Thus all the required quantities are known in the group height 
integral of Equation (12) . Although the integration cannot in 
general be performed analytically, it is nevertheless relatively 
simple, when N(h) is known, to compute p'(f) by a numerical inte- 
gration technique. The basic problem involved in the analysis of 
an ionogram is to perform the opposite conversion, namely to derive 
the N(h) function from a knowledge of the p'(f) function. It is, 
however, not generally possible to invert analytically the group 
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height integral. The method used is to find a general model for 
N(h), with many adjustable parameters, which will satisfy Equation 
(12) for selected values of the p'(f) function. The number of 
parameters which can be determined is the same as the number of 
p ' (f) values selected for the analysis. With the lamination 
technique used in this report the N(h) function is represented by 
a number of points (Nj,hj) connected by simple analytic curves. 

More specifically, in a given height interval (hj to hj), the 
profile is assumed to be of the form: 

h = h j;L + Fj(N) (13) 

where F.(N) is a simple analytic function of density and where 
Nj would be the density at h j . For the parabolic-in-log N method, 
the electron density profile is assumed to consist of k laminations 
as shown in Fig. 4. The first lamination (h^, h^) which begins at 
the satellite height h^ and extends down to the height h^, is 
assumed to be linear-in- log N. All other laminations, such as 
(hj hj) between heights hj ^ and hj are assumed to be para- 
bolic-in-log N with continuous slopes at the boundaries. A de- 
tailed treatment of the extraordinary trace analysis based upon 
parabolic-in-log (N) laminations is given in Appendix C. The 
lamination method, however, is more readily visualized in terms 
of the following simplified set of conditions. 

Assume that 1) the p'(f) function under analysis is the 
ordinary trace of a topside ionogram; 2) the trace is defined for 
fN < f < fm, where fN is the plasma frequency at the satellite 

and fm is the maximum frequency at which (0) reflections are 
obtained; and 3) the corresponding N(h) profile will be approxi- 
mated by a succession of linear segments (i.e., within a given 
lamination h varies linearly with N) . The procedure which will be 
described is such that the laminations are derived from a sequence 
of (p f , f) values selected prior to analysis. In the sequence 
(p’j, f j) , f increases monotonically with j and (Pq, f Q ) corresponds 
to zero range on the (0) trace, i.e., p^ = 0 and f^ = fN. 


- 10 - 



For the assumed linear lamination, Equation (13) becomes: 


h = h j _ 1 + aj (N-N j _ 1 ) (14) 

from which 

dh = a. dN (15) 

J 


Writing Equation (12) in terms of the laminations shown in 
Equation (14) and changing from the variable h to the variable N 
yields : 


p'(fj) 


i=j 


i=l 


!»Nj 


n' dN 


N 


i-1 


(16) 


where N. is the density at which reflection occurs at the fre- 

quency f i.e., N. is given by Equation (10). The right hand 
J J 

of Equation (16) represents an integration over (j) laminations. 
For the first lamination, i.e., the lamination nearest to the top- 
side sounder: 

P'(fi) = a x r"l n'CN.fj.B.e) dN (17) 

M 0 

The integral in Equation (17) can be evaluated with adequate 
accuracy by assuming B and 9 to be constants and equal to their 
values at the satellite. This approximation can be made because 
the ordinary ray group height integral is not very sensitive to 
the small variation in B and 0 which occurs within the altitude 
range of a typical lamination. With this assumption, a^ is com- 
pletely defined by Equation (17) and from Equation (14) it follows 
that 


h 1 = h Q + a x (Nj-Nq) 
where 

11 q = height of the satellite 

Nq = electron density at satellite (given by Equation (10) 
and letting f =* fN at the satellite). 
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For the second lamination: 


p’(f 2 ) = a L r N l n’(N,f 2 ,B o ,0 o )dN + a 2 T N 2 n ' (N,f 2 ,B 1 , B^dn (18) 
N~ N , 


where N 2 is related to f 2 by formula (10) . It should be noted 
that the integral associated with a^ is now for the frequency 
f 2 , and also that the value of B and 0 used in the second integral 
corresponds to the altitude h^ (which is obviously a more correct 
estimate of B and 9 in the second lamination). Equation (18) yields 
a 2 since this is the only unknown quantity, and consequently: 

h 2 = h l + a 2 ^ N 2~ N 1^ 

This step-by-step procedure is continued until the entire profile 
is determined. The relatively simple procedure described above 
has been used for the analysis of ionograms obtained from ground- 
based sounders (Jackson 1956). An additional simplification which 
is permissible in the (0) ray analysis of ground-based soundings 
is to use constant values of B and 9 at a given location, since 
the ionograms correspond to a relatively small altitude range 
(typically 100 to 300 km). In topside N(h) analysis, it is usually 
also permissible to treat 9 as a constant (value of 9 at the 
satellite for a given ionogram) , but the altitude variation of B 
must be taken into consideration. The techniques used for the re- 
duction of topside ionograms are basically refinements of the 
simple lamination concept outlined above. Section IV gives a 
general discussion of these refinements, and Section V indicates 
the improvements in accuracy resulting from these refinements. 

A more detailed discussion of the N(h) reduction technique is 
given in Appendices A, B, and C. 


IV. DISCUSSION OF THE LAMINATION METHOD 
Selection of Data Points 

From the preceding discussion of the lamination procedure, 
it was seen that the number of laminations obtained is determined 
by the number of p'-f values selected for the ionogram reduction 
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The choice of p'-f values is not particularly critical, provided 
these data points are sufficient in number and adequately distributed. 
Methods developed for the selection of data points take into con- 
sideration, not only the computer time required, but also the 
available scaling facilities. For example, at the Goddard Space 
Flight Center (GSFC) , ionograms are scaled with a graphic digitizer, 
which yields several hundred p'-f values per ionogram. Since this 
number of data points is one order of magnitude greater than required 
for N(h) analysis, the automatic analog-to-digital conversion is 
followed by a selection process. The selection process is a computer 
operation, done according to a fixed rule (GSFC criterion) and 
designed to provide enough data points in all cases (typically 
25 points), while avoiding excessive computer time. In many cases 
the GSFC criterion (Jackson 1967, p. 6 of appendix B) yields more 
points than necessary, but this slight disadvantage is more than 
compensated by the simplicity of the system. The GSFC criterion 
was used to scale the ionograms discussed in this report. 

Lamination Model 

Closely related to the number of points used in the analysis 
is the model used in the lamination technique. If a large number 
of points is used, then the resulting profile will be fairly ac- 
curate, even with the simple linear model used in the earlier example. 
However, the number of numerical integrations required increases as 
the square of the number of p'-f values used in the calculations. 

With a more elaborate model fewer points can be used in the cal- 
culation with an attending reduction in the computer time required, 
or a greater accuracy can be achieved using the same number of p'-f 
points in the calculations. Actually, in the topside ionosphere 
the electron density profile is represented more accurately by i 
succession of exponential segments, (Fitzenreiter and Blumle, 1964), 
i.e., the height increments are almost linear in log N, namely: 

Ah = a. (In N - In N._,) (19) 

J J 
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One objection to the linear-in-N or linea r-in-log N representation 
is that the assumed profile has discontinuous derivatives at each of 
the lamination boundaries. This difficulty is readily overcome by 
assuming that the height increments are parabolic with continuous 
slopes at the boundaries. The parabolic-in-log N assumption 
(Paul and Wright, 1963: Doupnik and Schmerling, 1965) assumes that: 

Ah = a. In N/N j _ 1 + b. [In N/N^] 2 (20) 

with the slope continuity yielding: 

Vi ■ + 2b i ln ‘WP (21) 

The number of laminations is determined by the number of scaled 
p’-f points, and the actual values of each p'-f pair determine the 
end points of the N-h laminations. However, since the analysis 
yields also the equation of each lamination, it is possible to 
calculate intermediate N(h) values to the same accuracy as that 
of the end points. Thus the analysis can yield either densities 
at fixed heights (such as at multiples of 100 km) or heights for 
pre-selected values of fixed densities. 

Evaluation of the Group Height Integral 

The following comments are concerned with the method used 
for evaluating an integral such as the one shown in Equation (17). 
There are three types of problems to be considered, the first one 
involves the parameters used in the integrand, the second one is 
concerned with the limits of integration, and the third one is the 
integration technique itself. 

In the example given earlier (0-trace, linear laminations) 
it was assumed that B was constant in each lamination. This 
yields a fairly accurate value of the height h^ and consequently 
of the altitude interval over which the integration is performed. 
Having determined the parameter a^ the altitude and the value of 
B are known for each value of N used in the numerical integration. 
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Hence the integral can be evaluated again, this time associating 
a more accurate value of B with each value of N used in the inte- 
gration. This will yield a slightly different value of a-^ and h^. 

The process converges very rapidly and after a couple of iterations 
there are no further significant changes in the final answer. 

The principle of iteration is slightly more complicated when 
extraordinary data are used in the calculations, since in this 
case it is not only the integrand which is affected but also the 
upper integration limit. If Equation (17) referred to a virtual 
height for the extraordinary ray, then the upper limit of inte- 
gration would be given by Equation (10) for the extraordinary 
trace (f = f ), namely: 

-L X 

N x = 12,400 f 1 (f 1 -2.8B 1 ) 

The value of B^ is not known, and it would have to be initially 
estimated by letting B^ = B^. Solving Equation (17) with this 
assumption would yield a fairly good estimate of h^ and hence B-^. 

The procedure could then be repeated, using not only a more ac- 
curate value of the integration limit N^, but also more representative 
values of B within the integrand. It is natural to anticipate 
that iteration should be more important for the extraordinary ray 
than for the ordinary ray. Further discussion of the iteration 
technique is given in Appendix C. 

The final poi.it is concerned with the integration technique. 

The problem which arises here is the fact that the integrand is 
infinite at the reflection point. Although it has been known for 
at least 18 years (Poeverlein, 1951; Shinn, 1951; Jackson, 1956) 
that this infinity can be removed by means of a suitable change of 
variable, the importance of this transformation was not fully ap- 
preciated until quite recently. Some of the early topside N(h) 
programs were based upon the assumption that a 16-point Gaussian 



integration technique could yield an accurate answer for an in- 
tegral such as the one appearing in Equation (17) . It turns out 
that this technique is both inefficient and inaccurate. Calculations 
(see Appendix B) for the (X) wave, longitudinal propagation, and an 
exponential lamination terminating at the reflection point, show 
that the numerical integration will be in error by 5.5% using a 
7-point Gaussian and by 2.5% using a 16-point Gaussian, whereas 
the error is less than 0.005% if the same integral is evaluated 
numerically after making the suggested change of variable and 
using only a 3-point Gaussian integration technique. Similar 
results are obtained when the propagation is not longitudinal. 

The accuracy of the (0) wave group retardation calculations 
can alsc be considerably improved by making a change of variables 
which keeps the integrand finite. A 3-point Gaussian is not ade- 
quate, however, except for low dip angles. As the dip angle is 
increased from 50 to 85 degrees, the number of points used in the 
Gaussian integration has to be increased from 3 to 16 in order to 
maintain an accuracy of better than 0.1%. Table 1 illustrates 
the errors for typical Gaussian integration and for various values 
of the magnetic dip angle. The errors shown in Table 1 are for 
the first lamination of the high density Bauer profile discussed 
in Section V. The density ratio N^/l^ for this lamination is 
1.05. Using the 3-point Gaussian, the errors are typically three 
times larger for the 10th lamination (density ratio = 1.18) and 
five times larger for the 27th lamination (density ratio = 1 26). 
Comparable results were obtained with the low density Bauer profile 
discussed in Section V. Similar calculations (for the bottom- 
side ionosphere) by Shinn and Whale (1951) led to the conclusion 
that their results became inaccurate when the dip angle was equal 
to or greater than 81 degrees. Becker (1960) did not give data 
beyond 80 degree dips in his tables. Present computer accuracies 
and integration techniques yield satisfactory resuxts up to at 
least 87 degrees. 
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DIP 

3-POINT 

7-POTNT 

12-Pt_.rr 

16-POINT 

0-50° 

— 

— — 





0 

o 

t> 

0.1% 

0.01% 

— 

— 

80° 

1.1% 

0.1% 

— 

— 

85° 

5 . 5% 

0 . 5% 

0.1% 

— 

0 

i> 

00 

17.5% 

1.5% 

0.3% 

0.1% (est.) 

89° 

200% 

10% 

1.5% 

0.5% (est.) 


Table 1 - Error in O-ray group height integrals for a 

lamination terminating at the reflection point. 
Blanks in the above table indicate that errors 
are much less than 0.1 percent. 

The integrand is kept finite at the reflection point by means 
of the following change of variables: 

2 

t = 1-X for the ordinary ray 

2 

t = i-X/(l-Y D ) for the extraordinary ray 
K 

where : 

X = 80.6 N/f 2 

Y^ = value of Y at reflection point 
Y = f H/f 

In the above commonly used expression for X, N is in electrons/cc 
and f is in kHz. 

The change of variable eliminates the infinity at the re- 
flection point for the following reason. The formulas giving 
n' near a giver reflection point (Jackson 1967, Appendix A, 
p. A7 and A8) are of the form: 

n' -- K/vT - kN 
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where K and k are constants appropriate for the propagation 

2 

mode used. By letting t = 1-kN, a typical group height in- 
tegral such as the one shown in Equation (16) becomes: 


1 n’dN 
i- 


2 /» t i 


* 1-1 


n’t dt 


At the reflection point, the integrand n’t is finite since: 


n’t = n v 1-kN * K 

V. DISCUSSION OF p'-f TO N(h) CONVERSION ERRORS 

Accurate analysis of topside ionograms require special care 
in a number of areas, such as the selection of data points, the 
choice of lamination technique, the use of iteration in the cal- 
culations, and the method used for the numerical evaluation of the 
group height integral. The importance of these considerations can 
be tested with ionograms corresponding to a known electron density 
distribution. The first test is to determine how well the top- 
side p'-f reduction technique (with all its present refinements) 
will reproduce the known N(h) profile. In subsequent tests various 

refinements are omitted and the resulting errors determined. 

xn 

The reference N(h) profiles (Fig. 5) used/ the error study 
were two theoretical topside N(h) distributions based upon a 
ternary ion mixture in diffusive equilibrium (Bauer, 1962). It 
is assumed that the ionograms were obtained from an altitude of 
3000 km (which corresponds to the apogee of Alouette II) and that 
the local electron densities were respectively 1000 el/cc for the 
low density profile and 5000 el/cc for the high density profile. 

For this theoretical situation it was assumed that the gyro-fre- 
quency at the satellite was 0.38 MHz and that the magnetic field 
varied with altitude according to an inverse cube law. The 
corresponding ionograms (X-traces) are shown in Fig. 6 and 7. 

To compute the ionograms, the profiles were divided into 
linear-in-log N laminations 5 km thick. The scaling criterion 
mentioned earlier was applied to the set of 520 p’-f values 


- 18 - 


thus obtained, yielding the points (open and solid) shown on 
the ionograms. The solid points illustrate a less detailed 
scaling in which approximately half of the data points would 
be eliminated. Based upon additional p'-f calculations, using 
10 km and 20 km laminations, it was concluded that the errors in 
the p'-f data (obtained with the 5 km laminations) were less than 
1 km. 

The normal GSFC procedure for reducing the (X) trace is to use 

the data points given by the scaling criterion, to assume laminations 

parabolic-in-log N, to iterate the height calculation for each new 

lamination until the results converge to within 0,01 km, and to 

integrate with a 3-point Gaussian after making the change of 
2 

variable t = 1-X/(1 -Yjj) . Using this procedure, the maximum error 
in altitude is about 1 km for both the high-density and the low- 
density profiles. This error is small compared to the errors due 
to uncertainties in the scaling of ionograms. Even on excellent 
ionograms the scaling error on the virtual heights is at least 5 km. 
Hence the recommended calculation procedures will not contribute 
significantly to the total error in the ionogram reduction process. 
Furthermore, even if scaling errors could be made negligible, a 
maximum error of 1 km would not be significant, particularly on 
a profile extending from 400 to 3000 km. Deviations from the 
recommended procedure, however, can sometimes introduce large 
errors as will be shown in the following discussion. 

1 , Selection of data points 

Satisfactory accuracy in the N(h) reduction of topside iono- 
grams can in general be achieved using 10 to 20 well selected p'-f 
values. For reasons discussed earlier, the GSFC criterion will in 
general yield a number of data points slightly greater than necessary. 
More importantly the criterion will very seldom yield an insufficient 
number of points. For example, on the low density profile the crit- 
erion yields 24 data points and a maximum error in altitude of 1 km. 

If the calculations are done with only 13 data points the maximum 
error becomes equal to 6 km, which is still relatively small. The 

error for these two cases is shown as a function of altitude on 
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Fig, 8 and 9 (graphs labelled: parabolic). These results indicate 
that the number of points given by the scaling criterion is not 
marginal, since it is possible to eliminate approximately half of 
these points and still achieve a satisfactory accuracy. 

2 o Choice of lamination model 

Performing the same calculations with the low density profile, 
but assuming laminations linear in log N, yields a maximum error 
of 20 km for the 24 point analysis, and 50 km for the 13 point 
analysis. The error for the linear lamination calculation is also 
shown in Fig. 8 and 9 (graph labelled: linear). Similar results 
were obtained for the high density profile. In this case the linear- 
in-log N lamination yielded a maxinmm error in altitude of 15 km 
using the 28 points provided by the scaling criterion, and the max- 
imum error was 40 km when only 15 points were used in the calculations. 

This leads to the conclusion that the parabolic lamination 
method yields results about 10 times more accurate than those 
obtained with the linear lamination method. Furthermore, the use 
of the linear lamination causes the profile to be too high. 

3 . Importance of iteration 

From the discussion in Section IV, it is seen that the con- 
cept of iteration applies only to the last lamination, i.e., the 
lamination nearest to the reflection point. Iteration is important 
only when the evaluation of this lamination is very sensitive to 
the height variation of the magnetic field. Thus iteration is 
particularly important for the analysis of the (X) trace when the 
corresponding electron densities are very low. Errors which can 
arise in the absence of iteration are shown in Fig. 10 for the 
theoretical ionogram of Fig. 6 (based upon the high density profile) 
and for two lamination models. It is seen that the largest errors 
occur when the parabolic-in-log N laminations are used. The 
error curve for the parabolic lamination also exhibits large os- 
cillations since an error in one lamination provides both an 
incorrect starting point and an incorrect initial slope for the 
next lamination. Repeating the above calculations with the ionogram 
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for the low electron density profile, leads to errors so large 
that the resulting profile is no longer monotonic. 

The parabolic-in-log N lamination leads to much greater 
accuracy, however, when iteration is performed (see Fig. 8 and 9). 
The linear-in-log N lamination provides a much better initial 
estimate of the lamination thickness and this estimate is used as 
the starting point for the parabolic technique. This procedure 
speeds up the convergence process (the 0.01 km convergence criterion 
being usually satisfied after 3 iterations) and also helps in pre- 
senting the parabolic calculation from going astray. 

4 . Choice of integration technique 

To show the importance of the change of variable in per- 
forming the integration, the ionogram for the low density profile 
was analyzed using the parabolic-in-log N lamination technique, 
iterating until Ah was less than 0.01 km, but omitting the change 
of variable. The integration was performed with both a 7-point 
Gaussian and a 16-point Gaussian. The error is shown as a function 
of altitude in Fig. 11. Even with a 16-point Gaussian an error of 
22 km can take place. The errors are such as to make the profile 
appear too low. 

5 . Comparable observations using actual Alouette II ionograms 

An error analysis was made on two Alouette II ionograms, 

taken from altitudes of 958 and 2873 km, with local electron den- 

4 3 

sities of 2.3 x 10 and 1.3 x 10 respectively. The analysis used 
a refined magnetic field model (Daniels and Cain, September 1965 
model), and it was assumed arbitrarily that the parabolic-in- log N 
techniques would yield a correct answer. The results, shown in 
Fig 12, indicate that the errors found with the high altitude 
ionogram were comparable to the errors found with the low density 
Bauer profile. Curve C for the high altitude case (Fig. 12) was 
smoothed; prior to smoothing it was similar to the corresponding 
curve on Fig. 10. Although a 7-point Gaussian integration was 
used to calculate curves C, the same results would be obtained 
with a 3-point Gaussian integration. It is also seen "<hat the 
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errors are much smaller for the low altitude ionogram. In part- 
icular, iteration is not nearly as important for a low altitude 
ionogram as for a high altitude ionogram. 

Errors obtained when an inverse cube law is used (based upon 
the correct value of B at h g ) have been investigated with these 
two test ionograms. The maximum error resulting from the use of 
an inverse cube field model was 10.4 km for the high altitude pro- 
files and 0.22 km for the low altitude profile. It was also found 
that no significant error is introduced if the magnetic variation 
is assumed to follow the inverse cube law within each lamination . 

6. Additional Considerations 

The preceding error study was limited to an investigation of 
errors arising from the numerical inversion of the group height 
integral. Systematic errors usually less than 30 km can also be 
present in the original p'-f data (see Franklin et al in this 
issue and "Comparisons between Topside and Ground-Based Soundings", 
Jackson, in this issue). The validity of the assumption of vertical 
propagation was not discussed in the present paper. It is shown 
in a separate paper ("Comparisons between Topside and Ground-Based 
Soundings", Jackson, in this issue) that the assumption of vertical 
propagation is a fairly good approximation in many cases, particularly 
for low altitude ionograms. 

Small errors (usually less than 5 kilometers) can also arise 
if allowance is not made for the fact that the sounder altitude can 
change by a few kilometers during a sounding. 

VI. CONCLUSIO N 

The techniques presently available for the analysis of topside 
ionograms can yield electron density profiles with an altitude 
accuracy of the order of a few kilometers, provided the p'-f data 
is free of scaling and of systematic errors, and provided the 
soundings correspond to vertical propagation. In actual practice 
errors in N(h) analysis are due primarily to the accuracy of the 
P'-f data and to deviations from vertical propagation. 
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APPENDIX A 


A- 1 


FORMULAS FOR GROUP AND REFRACTIVE INDICES 


1 . Basic Formulas 

The group refractive index n* for a radio wave of frequency 
f is defined as the free space velocity of light divided by the 
group velocity of the wave. The fundamental formula for n’ is: 


n’ 


f on 

Jf 


(A-l) 


where n is the real part of the refractive index of the medium. 
The index n is given by the well known Appleton Hartree formula, 
namely: 


n = 1 - 


1 - 


where 


l " 2 

T 

2(1-X) 
2 


V 


Y 2 
2(1-X) 


+ Y, 


X = N/(12,400f 2 ) = 


(A-2) 


N = electrons/cc 

fN = electron plasma frequency 

f = frequency (MHz) 

Y = fH/f 
fH = 2 . 8B MHz 

B = (terrestrial) magnetic induction in gausses 
Y t = Y Sin 0 

Y L = Y Cos 0 

$ = angle between propagation vector and magnetic field 

+ = positive sign in front of square root is for ordinary ray 
negative sign is for extraordinary ray. 
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A-2 


For the vertical propagation assumed in the analysis: 

Y T = Y Cos e 
Y l = Y Sin 6 

where 

c = angle of magnetic dip 

The reflection conditions are X=1 for the ordinary ray (except 
when 9 is exactly 90°) and X=l-Y for the extraordinary ray. 

The evaluation of i S fairly complicated since n is a 

function of X and Y and both of these parameters are functions 
of f. The calculations are simplified considerably, however, 
when 9 is either 0 or 90 degrees, i.e., when the earth's magnetic 
field is either perpendicular to (transverse propagation) or paral 
lei to (longitudinal propagation) the direction of the wave pro- 
pagation. For these two limiting cases it is also possible to 
evaluate analytically the group height integral fn dh and hence 
check the accuracy of the numerical integration technique which 
is needed for the general case (i.e., 9 neither 0, nor 9C degrees) 
2 . Special Cases 

a. Transverse Propagation (6=0; Y T =Y ; Y L =0) 

In this case formulas (A-2) and (A-l) give: 
for the ordinary ray: 

n = J 1-X (A-3) 

n’ = 1/n (A-4) 

for the extraordinary ray: 

n = / 1 - (A-5) 

1-X-Y 2 
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A-3 


n ' 


1 + 


XY‘ 


(1-X-Y 2 ) 2 

b. Longitudinal Propagation (6=90°, Y L =Y, Y^,=0) 
In this case formulas (A-2) and (A-3) give: 
for the ordinary ray: 


(A-6) 


n = jl - X/( 1+Y) 


r _ ± 

n 


XY 


1 - 


2(l+Y) J 


for the extraordinary ray: 


n = /I - X/(l-Y) 


? _ _ 


n ' = 


1 + 


XY 


2 ( l-Y; J 


(A-7) 

(A-8) 

( A - 9 ) 

( A- 10 ) 


3 . Doupnik's Formulas for the Group Index . 

The group index can be evaluated by substituting into Eq. 

(A-l) the value of n given by Eq. (A-2) . The resulting formula 
for the group index is complicated and it does not lend itself 
to accurate numerical calculations. An expression for n’ , which 
is compact and also well suited for numerical analysis, can be ob- 
tained by making a substitution suggested by J. R. Doupnik (private 
communication), namely: 
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A-4 


tan O' 


2Y t ( 1-X) 


(A -11 ) 


and evaluating equation A-l yields: 




(A - 12) 


where 




(A - 13) 


Coso 


S = 1+Y, -i — for the ordinary rav 

L 1+Sino J 


e = -1 


for the ordinary ray 


S = l-Y 


l+Sitio 
L Cos o 


for the extraordinary ray 


e = + 1 for the extraordinary ray 

The above formulas are not valid for Y^ = 0. In this case form;] 
(A-3), (A-4), ( A —5 ) and (A-6) are used in lieu of (A-ll). (A-12 ) 
and (A-13), 
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APPENDIX B 


B-l 


ACCURACY OF GROUP HEIGHT INTEGRAL CALCULATIONS 

h 2 

The evaluation of P n'dh requires both accurate values 

hi 

of the group index and an accurate integration technique. The 
accuracy of the n’ calculations can be checked against Becker's 
tables (1960) . The integration technique is critical only for 
the last lamination, i.e. , for the lamination which includes 
the reflection point (where n' becomes infinite). Hence the 
method of integration must be checked with an integral of the 
form: 

I - f hR n’dh (B-l ) 

h R-Ah 


where the upper limit of integration h^ is the reflection point. 
The accuracy of the numerical integration techu aue can be invest- 
igated in two different ways. One way is to perform the numerical 
integration several times, each time increasing the number of 
sampling points L. If the results become constant for L>K, then 
K sampling points yield an accurate answer. Another method is 
to examine one of the special cases when the integration can be 
performed analytically. The test discussed here v* s done using 
the n' function for the extraordinary ray, longitudinal propa- 
gation and constant Y. This yields an integral which can be 
evaluated analytically. A linear in log N lamination was used 
(see Eq. 19) giving: 


and 


h = hj_^ + (In N - In Nj j) 
,, dN dX 

dh ~ a j U " a j T 
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Changing to the variable X in Eq. ( B-l) gives: 

I = a . f Y §dX 

Since is outside of the integral, it can be assumed to be 
unity for the purpose of checking the integration technique. 
In this case n* is given by Equation (A-10) , which can be 
written 


n’ 


1+bX 

yi+aX 


where a = » and b = tr 

i_Y 2(1-Y) ^ 

For any specified value of Y, the quantities a and b are 
constants. Substituting n’ from Er . (B-3) into Eq. (B-2) 
yields : 


iiSL_ | dx 

X/T+aX ' 

The recommended substitution for the numerical evaluation of 
Eq. (B-4) is to let: 



1-X/ (1-Y) 


giving: 


2tdt = - dX/ (1-Y) 


The integration limits are; 

t - (l _ _ / c 

t l " U 1-Y ' " «/ 1-Y 


and 


t 2 0 


(B-2) 


( B-3) 


(B-4) 
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B-3 


After making the above change in variable, and noting that for 
this special longitudinal case t = \l+aX, Eq. (B-4) becomes: 


I 




.v T-y 


2 (1-Y) (• 


1+bX 

X 


) dt 


(B-5, 


where: X — (l-t^) (1-Y) . 

The analytical solution of Eq. (B-4) is: 

X=l-Y 
X=l-Y-e 


T . 2Wl+aX , ,1-k/ 1+aXv 

I = | — In ( ) 

a 1-v 1+aX 


In the above expression, the ratio of the integration limits 
(l-Yj/d-Y-e) is the ratio of the values of X at the two boundaries. 
Since X is proportional to N, this ratio is also the electron 
density ratio Ng/N^, where is the initial (or minimum) density 
and is the last (or maximum) density for the lamination. The 
numerical and analytical integrations were carried out for valr.es 
of Y ranging from 0.1 to 0.9 (for routine analysis Y will seldom 
exceed a value of 0.90) and for electron density ratios (l-Y)/(l-Y-e ) 
ranging from 1.01 to 2=0. The numerical integration was performed 
with the indicated change of variable using a 3 point Gaussian 
technique, and also without making a change of variable with a 
7 point and a 16 point Gaussian technique. 

The errors arising from the various numerical integration 
techniques are shown in table B-l. The error is roughly constant 
and quite substantial for each Gaussian integration where no change 
in variable was made (4.7 to 5.8 percent with the 7 point Gaussian 
and 2.1 to 2.6 percent with the 16 point Gaussian). When the 
change in variable is made, the error is less than 0.1 percent 
for all conditions considered and only 3 points are required in 
the Gaussian integration. This remarkable improvement in accuracy 
is obtained because the integrand in Eq. (B-5) is a very slowly 
changing function. For the special case investigated here (9=0) 
the curves representing the integrand as a function of t are 
parallel to each other as the parameter Y is changed. This is 
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B-4 


readily seen for the case 9=0; since then: 

2 (1-Y) n ’ t/X = — \ + JL, 

The integrand obtained for the general case (9^0) has a similar 
slow variation as a function of t as can be seen from Fig. B-l, 
and consequently the data shown in table B-l is also indicative 
of the improvement in accuracy which can be obtained for arbitrary 
values of 0 when a change in variable is made 0 

In order to apply these results to the actual N-h analysis, 
it is necessary to calculate the density ratios which are obtained 
with the GSFC criterion for the selection of data points. The 
density ratios depend upon fH and fx at the satellite and also 
upon the frequency increments specified by the criterion. Typical 
ratios are shown in Fig. B-2. It is seen that the ratios are 
always less than 1.6. Hence the maximum error obtained using 
the 3 point Gaussian with change of variable is less than 0.020. 
Thus for the extraordinary ray this method is at least two (and 
typically three to four) orders of magnitude more accurate than 
the 16 point Gaussian with no change in variable. 
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B-5 



7 

POINT GAUSS' 

::AN - 

- NO CHANGE 

IN VARIABLE 


Y 

1.01 

Density 

1.10 

Ratio ir. Lamination 
1.40 1.60 

2.00 

.1 

5.79 

5.63 


5.23 

5.02 

4.70 

.3 

5.79 

5.66 


5.29 

5.11 

4.82 

.5 

5.79 

5.68 


5.39 

5.23 

4.98 

.7 

5.79 

5.72 


5.51 

5.39 

5.21 

.9 

5.80 

5.77 


5.69 

5.64 

5.56 



16 

POINT GAUSSIAN 

- NO CHANGE 

IN VARIABLE 


.1 

2.63 

2.56 

2.38 

2.28 

2.14 

.3 

2.63 

2.57 

2.41 

2.32 

2.19 

.5 

2.63 

2.58 

2.45 

2.38 

2.27 

.7 

2.63 

2.60 

2.50 

2.45 

2.37 

.9 

2.64 

2.62 

2.50 

2.56 

2.53 



3 POINT GAUSSIAN 

- CHANGE 

IN VARIABLE 


.1 

0.000025 

0.000044 

0.0041 

0.014 

0.065 

.3 

0.000025 

0.000038 

0.0036 

0.013 

0.058 

.5 

0.000025 

0.000031 

0.0029 

0.011 

0.048 

.7 

0.000025 

0.000021 

0.0021 

0.008 

0.035 

.9 

0.000025 

0 „ 000007 

0.0008 

0.003 

0.015 


Table B-l. Percent error in the evaluation of the group 

!-Y n T 

height integral — dx as a function 

of Y and of the electron density ratio 
( l-Y)/(l-Y-e) for various integration techniques. 
Calculations were performed for the case of 
longitudinal propagation. 
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APPENDIX C 


N (h) CALCULATIONS USING THE EXTRAORDINARY TRACE AND PARABOLIC- 
IN LOG (N) LAMINATIONS. 


1. General Procedure 

It was indicated earlier (Section II of the text) that most 
of the topside N(h) profiles have been derived from the analysis 
of the extraordinary ray echoes. The principles of the analysis 
technique (lamination method) were introduced in Section III and 
discussed further in Section IV. In Section V it was shown that 
very accurate results can be obtained if the laminations are 
assumed to be parabolic-in-log (N) with continuous slopes at 
the lamination boundaries (see Fig. 4). The present appendix 
is devoted to a detailed description of the (X) data analysis 
using parabolic-in-log (N) laminations and assuming vertical 
propagation. The notation used is the same as in Section III, 
and for convenience a few formulas previously given are repeated 
here . 

*t hi 

The j— - lamination is defined by: 


where 


j-l + a j ln + b j< ln 

(C-l) 

N._, 

a . , + 2b . , ln 
j-l j-l Nj^ 

(C-2) 


The above formulas apply to all laminations except the first 
one, which is assumed to be linear-in-log (N) , i.e. , b^ = 0. 

The calculations also require that all expressions in (C-l) and 
(C-2) be finite. This requirement will be met only if the N(h) 
function is monotonic. 

In the N(h) calculations, the laminations are determined 

one at a time, in the order (h^, h^) , (h^ , hg) , (hg, h^) etc. 

Hence all laminations between the altitudes h_ and h. , are known, 

th 0 3 ~ 1 

when the calculation of the j — lamination is performed. Let f . 

represent the frequency of the extraordinary wave reflected from 

the unknown height h.. The virtual height h’ . is given by: 

3 J 


h h 

r o , i 

p f . = j n ' dh + 

J " h l 


n ' dh + 


+ 


h . 
J-l 


n 'dh 


(C-3) 
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h . 
J 


C-2 


DP + I J 1 n’dh 

• h . 

3 


(C-4) 


where DP represents the delay in the previous laminations. The 
calculation of DP presents no special problem, since it involves 
calculation of integrals in which all parameters are known. The 
process for calculating DP will be illustrated later. For the 

present discussion (calculation of j— lamination) it will be 
assumed that DP is known. 

2 „ Initial Calculation of .j— lamination (Constant Y) 

The reflection conditions at the bottom of the j— lamina- 
tion can be written in terms of the X and Y notation of Appendix A 
as follows: 


X. = 1 - Y. 
J J 


(C-5) 


where 


(fH), 


Y. = 
J 


(C-6) 


f . 
J 


Equation C-l can be expressed in terms of X by noting that 


N 


N 


3-1 


x j-i 


( C — 7 ) 


Making the substituting indicated by Eq. (c-7) in Eq . (c-l) 
and differentiating Eq. (C-l) yields: 


dh = [a + 2b In («£—)] ~ 

J J j-i X 

Substituting the above value of dh in Eq. (C-4) yields: 


(c-s) 


P' . = DP + a . 
J J 


x j-i 


x j-i 


f 2— ax + 2b. [ 5- ln(^— ) dX fc-9) 

"X. x J J x. A A j-i 
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Letting t" = l-X/(l-Y.) and noting that dX = -2(l-Y )t dt 

3 J 

gives 

P’ j = DP - a. J tj_1 2(1-Y.) ^ t dt 

- b. f j" 1 4 ( 1-Y . ) Si ln(^— ) t dt ( C-10) 

,1 3 ^ A j - 1 


To simplify the notation of Eq. C-10, the parameter X was retained 
in the formula. In Eq, C-10, however, X is the following function 
of t : 


X = (1-t 2 ) (1-Yj ) 


(C-ll) 


Similarly, since n' is a function of X, Y and 6, the values of 
X entering into the calculation of n' are those given by Eq. C-ll. 
Since the value of Y. is not known initially, it is necessary 

J 

to use an estimated value of Y, to compute the integrals in Eq. 
(C-10). One method is to let Y. = Y._, anci to assume that Y is 
constant within the lamination and equal to Y. , . All the para- 

J 

meters are then known in Eq. (C-10) and the integrals can be 

evaluated. Representing the integrals associated with a. and 

J 

b. by S and S respectively, yields: 

J A D 

bj = (DP - p'j - aj S A )/S B (C-12) 


The above initial estimate of Y. can cause the iteration pro- 

J 

cess to fail. For this reason a different procedure is used as 
indicated in section 4 of Appendix C. 


3. Iteration With Variable Y 


The approximate answer thus obtained for the lamination 
can be refined by an iteration process, in which Y. is computed 

J 

for the calculated value of h. and in which the values of Y 

3 

entering into the n* calculations are computed by assuming that 

Y decreases from the value Y. , at the top of the lamination to 

J 

the value Y at t he bottom of the lamination according to: 

J 
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Y — — — 

K 

(C-13) 

(R 

+ h) 3 


The constant K and the 

■ earth's radius R in Eq. (C-13) 

can be 

eliminated by making use of the boundary conditions: 


Y.(R ♦ h/ = 

Y i-i < R + h i-i> 3 = K 

(C - 14 ) 

To express Y in terms 

of Yy Yj_ 1 and h, Eqs. (C-13) 

and (£-14 

are written: 



R + h 

k/y 1 / 3 

(C-15) 

R + h. 
J 

k/y . 1//3 

' 3 

(C- 16) 

R + h.^ * 

K/Y. ^ 1/3 

(C-17) 

from which 



h - j 

K( Y -1 / 3 - Yj' 1 ^ 3 ) 

(C-18) 

h. , - h. 
3-1 J 

K ( Y j_l 1//3 " Y j~ 1//3 ^ 

(C-19) 


The constant K is eliminated by dividing Eq. (C-18) by Eq. (C-19). 
Solving the resulting equation for Y yields: 


Y. 


Y /. (h-h.) .3 

1 ♦ r(y^-) 1 / 3 -i] (h . J.) } 

.1-1 ' .i-i .r 


(C-20) 


The values of h to he used in Eq. (c-20) should be those corres- 
ponding to the values of X. It is therefore important to note 

that the b. obtained from Eq. (C-12) was based upon the assumption: 
J 


X. 

J 


'-Vi 


(C-21) 


Hence to be consistent the altitude h. must be calculated according 

3 

to the formula : 
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h. = h. , + a. ln( Y ^- ) + b. [in 1‘ 

J J-l J Vi J X j-1 


(C -22) 


where X . = 1-Y . 

J J-l 

If Y. is now redefined in terms of h., (giving X. -> 1-Y.) 

J J J J 

Eq. (C-22 ) will no longer yield the same value of 

h. at the reflection point. In order to proceed with the 
J 

iteration (and in particular make sure that Eq . (C-2Q) will keep 

Y between Y. , and Y.), we must either redefine h. in terras of 
J-l J J 

X. = 1-Y., or recompute b. so that Eq . (C-22) gives the same 
J J J 

value of h. for the X.. In the GSFC program, the computed 

J J 

value of h . Is preserved. The parameters Y., X., and b. are 
j 1 j’ j’ j 

redefined prior to the iteration process. Iteration is continued 

until the successive values obtained for h^ agree to within 0.01 km. 

When this happens the difference between the b^ computed from 

Eq„ (C-12) and the b. reevaluated prior to iteration becomes 

J 

insignificant. When the desired convergence is achieved, the 

final value of b . is the value computed from Eq . (c-12) and the 
J 

final value of h. correspond to the Y. and X. obtained from the 
J J J 

previous calculation. Hence the final compromise is made on 

fHj which is actually computed (and stored in the program) for 

an altitude slightly different than h. (the altitude difference, 

. J 

however, being less than 10 meters). 

Returning to the postponed d iscussion of DP of Eq. (C-4), 

it is seen that DP involves a summation of integrals identical 

to those of the j — lamination, except that the limits are 

different, but known when the j— lamination is calculated. 

These integrals are also evaluated making use of the change of 

variable t 2 = 1 - X/(l-Y.), and using Eq. (C-20) to vary Y within 

J 

the lamination. 

For the first lamination, which is assumed linear-in-log N, 
the coefficient _b is zero and Eq. (C-10) is solved for the co- 
efficient a. Actually the first calculation made on the j — 
lamination is also based upon a linear-in-log N method; this 
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yields a good estimate of the value of hj , which is then used 
as the starting point for the iteration using the parabolic-in- 
log N technique. 

4. Convergence of the iteration process . (Initial value of Yj) 
The p'-f to N(h) conversion discussed in this paper requires 
that the p’-f trace be well defined and equivalent to a continuous 
curve. This requirement is easily met on topside ionograms, be- 
cause in the topside ionosphere N usually decreases monotonically 
with h. Such an N(h) distribution is implicitly assumed in the 
parabolic-i n-log N lamination technique and also in the trans- 
formation t^ = 1-X/(1-Y.). 

J 

The existence of a continuous trace on a topside ionogram 
usually implies that N decreases monotonically with h. Exceptions 
to this general rule are theoretically possible for the X trace. 
The reflection condition: 

fz - SU Vz 4 (fN) 2 + (fN) 2 


shows that echoes can be continuously received as a function of 
fx, provided the quantity S~4(fN)^ + fH^J increases monotonically 
with distance below the sounder. Since fH increases as a function 


of depth below the sounder, this increase could overcompensate 

a decrease in N (Lockwood, private communication). In the present 

N(h) reduction technique, such a situation is readily detected 

and it results in a rejection of the ionogram. 

Even if the profile decreases monotonically with altitude, 

the iteration process (based upon the initial assumption Y. = 

J 

Yj_^) may not converge if the profile is sufficiently steep. 

The permissible values of gyrof requencies at the bottom of the 

lamination range from the value at the top of the lamination 

f H . , to a maximum value fH.„ corresponding to the minimum 
j-i JM 

permissible increase in density. Thus, if we define the mir imum 
permissible increase in density as 0.1 percent, the quant it / 
fHj M is defined by: 


fx j (fx j - f V - (1 - 001) - fH j-i> 


(C-23) 
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True height calculations yield different answers for h. as 

J 

fHj is varied from fHj_j to fH . The basic purpose of the iteration 
process is to find from all the permissible values of hj , the par- 

J 

ticular height h- at which the f H . used in the true height calcula- 
J J 

tion is the same as the actual value of fH . at the altitude h.. The 

J J 

above iteration process is illustrated by Fig. C-l, which shows 

how true-height calculations are influenced by the assumed value 

of fH at the bottom of a lamination. The curves Lj , L 2 and L 3 

show the heights h^ at the bottom of the first lamination as a 

function of fH^ over the permissible range of fH^ for various 

assumed values of p^ on the extraordinary trace. The assumed 

conditions at the top of the lamination (i.e. , at the satellite 

altitude) were hS= 2000 km, fHs = 0.45 MHz and fxS = 0.50. At the 

bottom of the lamination, it was assumed that fi = 0.52 MHz. An 

inverse-cube variation was assumed for fH. For p T ^ = 1800 km, the 

resulting heights are given by the curve L^. Curve M shows the actual 
gyrofreque’ cy at the altitudes shown. The correct height is 
therefore given by the intersection of curves and M. If the 
initial assumption is fH^ = f HS , the first value of h.^ would 
correspond to A^. The iteration would then be performed using 
fH at Bj . which is the correct fH at the altitude of A^. This 
would yield a new value of h^, namely that corresponding to C^. 

The second iteration is performed using the fH value at 
It is seen graphically that the process converges to the point 
0^. The same process, however, will not converge on curve Lg , 
although there is a solution at Og . If the calculation is 
started with the maximum permissible value of fH, the initial 
height will be at on curve and at on curve L,. It is 
seen that the process will then always converge, if there is 
a solution. Curve L 3 is an example in which there is no solution 
such that the density at the bottom of the first lamination is 
greater than the density at the satellite. 
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In view of the above consideration, the GSFC program performs 
the initial calculation of the laminations using the value of 
fHj defined by Eq. (C-23) , provided fHj M is less than 1.15 fH^^. 
This upper limit was set for fHj M because at lower altitudes where 
the density increases rapidly with depth, Eq. (C-23) leads to un- 
reasonably high values of fHj M . 
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FIGURE CAPTIONS 
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Fig. 
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Fig. 

Fig. 

Fig. 

Fig. 

Fig. 


1. Typical Group Velocity Curves. 

2. Example of an Alouette II ionogram with excellent Z, 

0 and X traces. 

3 Electron densities obtained from the ionogram of Fig. 3, 

doing independent N(h) analysis on each of the Z, 0 and 
X traces. The points shown in the graph were selected 
from the computed points. For the sake of clarity, over- 
lapping points were omitted, except at the ends of the 
profile . 

4. Lamination Model. 

5. Ionospheric profiles used for error studies. 

6. Ionogram corresponding to the high density profile of 
Fig. 5. Points shown are those obtained using the 
author's scaling criterion. 

7. Ionogram corresponding to the low density profile of 
i-Lg. 5. Points shown are those obtained using the 
author’s scaling criterion. 

8. Relative accuracy of parabolic-in-log N and linear-in- 
log N techniques for the 24-point analysis. 

9. Relative accuracy of parabolic-in-log N and linear-in- 
log N techniques for the 13-point analysis. 

10. Errors introduced when iteration is not used. Graphs 
shown are for the linear-in-log N and for the parabolic- 
in-lcg N methods. Point A is common to both graphs 

and computed using the linear-in-log N method. The 
error on the next point is much smaller using the 
linear (B^) than using the parabolic-in-log N method (Bg) - 

11. Errors introduced when the group height integrals are 

o 

evaluated without making the change of variable t = 1 - 
X/(1-Y r ). 
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Figo 12. Errors introduced in ionogram analysis when either the 
integrand is infinite at the reflection point, or when 
iteration is not performed. The calculations were 
based upon laminations parabolic-in- log N. Hs is the 
satellite altitude. 

Fig. B-l. Behavior of integrand in extraordinary ray group height 

o 

integrals as a result of the change of variable t = 
1-X/U-Y) . 

Fig. B-2. Size of laminations (in terms of density ratios) re- 
sulting from the author's scaling criterion for various 
values of gyrof requencies and densities at the satellite. 
The gyrof requency was assumed to be independent of al- 
titude, which yields an overestimate of the ratios. 

Fig. C-l. Computed altitude for the bottom of a lamination as a 
function of gyrof requency . The minimum gyrof requency 
is the value at the satellite, and the maximum gyro- 
frequency is based upon an essentially constant density 
within the lamination. The actual gyrof requency as a 
function of altitude is also indicated. The correct 
value of gyrof requency at the bottom of the lamination 
is the intersection of the two curves. The purpose 
of the iteration process is to find this intersection 
point . 
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Fig. 5 
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Fig. 11 
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Fig. 12 



A: 16 PT. GAUSS, 3 ITERATIONS 
NO CHANGE IN VARIABLE 

B: 7 PT. GAUSS, 3 ITERATIONS 
NO CHANGE IN VARIABLE 

C: 7 PT. GAUSS, NO ITERATION 
CHANGE IN VARIABLE 
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